Symmetry energy of dilute warm nuclear matter 
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The symmetry energy of nuclear matter is a fundamental ingredient in the investigation of exotic 
nuclei, heavy-ion collisions and astrophysical phenomena. New data from heavy-ion collisions can 
be used to extract the free symmetry energy and the internal symmetry energy at subsaturation 
densities and temperatures below 10 MeV. Conventional theoretical calculations of the symmetry 
energy based on mean-field approaches fail to give the correct low-temperature, low-density limit 
that is governed by correlations, in particular by the appearance of bound states. A recently devel- 
oped quantum statistical (QS) approach that takes the formation of clusters into account predicts 
symmetry energies that are in very good agreement with the experimental data. A consistent de- 
scription of the symmetry energy is given that joins the correct low-density limit with quasiparticle 
approaches valid near the saturation density. 

PACS numbers: 21.65.Ef, 05.70.Cc, 25.70.-q, 26.60.Kp, 26.50.-^x 



The symmetry energy in the nuclear equation of state 
governs phenomena from the structure of exotic nuclei 
to astrophysical processes. The structure and the com- 
position of neutron stars depend crucially on the density 
dependence of the symmetry energy [l[. The symmetry 
energy characterizes the dependence of the nuclear bind- 
ing energy on the asymmetry S = (TV — Z) /A with Z and 
N the proton and neutron numbers, and A = N + Z. As 
a general representation of the symmetry energy coeffi- 
cient we use the definition 

Esym[n,T) = E{n,0,T), 

(1) 

where E(n, S, T) is the energy per nucleon of nuclear mat- 
ter with density n, asymmetry S, and temperature T. At 
low density the symmetry energy changes mainly because 
additional binding is gained in symmetric matter due to 
formation of clusters and pasta structures [2| • 

Our empirical knowledge of the symmetry energy near 
the saturation density, no , is based primarily on the bind- 
ing energies of nuclei. The Bethe-Weizsacker mass for- 
mula leads to values of about i?sym('T-o, 0) = 28 — 34 MeV 
for the symmetry energy at zero temperature and satura- 
tion density no ~ 0.16 fm~'^, if surface asymmetry effects 
are properly taken into account [1] . 

In contrast to the value of EgyminojO), the variation 
of the symmetry energy with density and temperature is 



intensely debated. Many theoretical investigations have 
been performed to estimate the behavior of the symme- 
try energy as a function of n and T. A recent review is 
given by Li et al. Q , see also 0, Q . Typically, quasipar- 
ticle approaches such as the Skyrme Hartree-Fock and 
relativistic mean field (RMF) models or Dirac-Brueckner 
Hartree-Fock (DBHF) calculations are used. In such 
calculations the symmetry energy tends to zero in the 
low-density limit for uniform matter. However, in accor- 
dance with the mass action law, cluster formation dom- 
inates the structure of low-density symmetric matter at 
low temperatures. Therefore, the symmetry energy in 
this low-temperature limit has to be equal to the binding 
energy per nucleon associated with the strong interac- 
tion of the most bound nuclear cluster. A single-nucleon 
quasiparticle approach cannot account for such struc- 
tures. The correct low-density limit can be recovered 
only if the formation of clusters is properly taken into 
account, as has previously been shown in Ref. in the 
context of a virial expansion valid at very low densities, 
and in Rcf. [1]. 

In this letter we employ a quantum statistical (QS) ap- 
proach which includes cluster correlations in the medium. 
It interpolates between the exact low-density limit and 
the very successful RMF approaches near the saturation 
density. We show that this picture is in agreement with 
recent experimental findings on £^sym at very low densi- 
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TABLE I: Temperatures, densities, free and internal symmetry energies for different values of the surface velocity as derived 
from heavy-ion collisions (cols. 2-6), from the QS approach (cols. 7-8) and selfconsistently with clusters (cols. 9-12), see text. 
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ties. 

Suitable approaches to account for cluster forma- 
tion are the nuclear statistical equilibrium (NSE) model 
0, cluster- virial expansions 0, and generalized Beth- 
Uhlenbeck approaches [loj . A thermal Green function 
approach that allows a generalization of the NSE model 
by introducing a quasiparticle description also for the 
bound states was already formulated some decades ago 



11| , but only recently analyzed with respect to the conse- 
quences for nuclear matter [l^l . In this QS approach the 
cluster correlations are described in a generalized Beth- 
Uhlenbeck expansion. The advantage of this method is 
that the medium modifications of the clusters at finite 
density are taken into account. In Ref. Q the thermo- 
dynamic properties of nuclear matter were derived using 
this approach. The formulation of Ref. Q is valid in the 
density and temperature range where the formation of 
light clusters with A < 4 dominates and heavier clus- 
ters are not yet important. The method requires a suf- 
ficiently accurate model for the quasiparticle properties, 
for which we ernploy a RMF model with density depen- 
dent couplings [l3| which gives a good description both 
of nuclear matter around normal density and of ground 
state properties of nuclei across the nuclear chart. In 
order to extend the applicability of this RMF model to 
very low densities, it has been generalized in Ref. Q to 
account also for cluster formation and dissolution. 

We note that at very low densities and temperatures 
below T ~ 1 MeV new phases may occur. In fact, the for- 
mation of a solid phase using Overhauser orbitals includ- 
ing a triple point 14 1 or Bose- Einstein condensation [lij 
have been suggested. However, in this letter we are con- 
cerned with experimental data which probe nuclear mat- 
ter at considerably higher temperatures. The low-T be- 
havior is an interesting issue for future studies. 

In the following we focus on finite temperatures and on 
the sub-saturation region n < uq. Experimental informa- 
tion is derived from heavy-ion collisions of charge asym- 
metric nuclei, where transient states of different density 
can be reached, depending on the incident energy and the 
centrality of the collision. In the Fermi energy domain 
symmetry energy effects have been investigated using ju- 



diciously chosen observables [4, 16- lij. 

Recently, the experimental determination of the sym- 
metry energy at very low densities produced in heavy ion 
cohisions of ^"Zn on ^^Mo and i^'^Au at 35 MeV per nu- 
cleon has been reported [2^. Results of this study are 
given in the first four columns of Tab. |T1 Note that as 
a result of an energy recalibration and reevaluation of 
the particle yields in different velocity bins these values 
are slightly different than those reported in Ref. [20j . 
The surface velocity Wsurf, i-e. the velocity before the fi- 
nal Coulomb acceleration, was used as a measure of the 
time when the particles leave the source under different 
conditions of density and temperature. Only values of 
Vsuri < 4.5 cm/ns are included here, since the system 
does not reach equilibrium for higher Wsurf, see Tab. I of 
Ref. The yields of the light clusters A < 4 were de- 
termined as a function of fsurf ■ Temperatures were deter- 
mined with the Albergo method [2l[ using a H-He ther- 
mometer based on the double yield ratio of deuterons, 
tritons, '^He and ''He, and are given in Tab. [J as the 
average for the two reactions. 

The free neutron yield is obtained from the free proton 
yield and the yield ratio of '^H to '^He. To determine 
the asymmetry parameter of the sources the total proton 
and neutron yields including those bound in clusters are 
used. The proton chemical potential is derived from the 
yield ratio of ^H to ^He. The corresponding free proton 
and free neutron densities are calculated, and the total 
nuclcon density is obtained by accounting also for the 
bound nucleons according to their respective yields poj . 
The total nucleon densities are of the order of 1/lOOth to 
l/20th of saturation density, as seen in Tab. HI 

An isoscaling analysis [221 been employed (as a 
function of Usurf ) to determine the free symmetry energy 
-Fsym via the expression a ~ 4i^sym A(Z/yl)^/r. Here a 
is the isoscaling coefficient determined from yield ratios 
of Z = 1 ejectiles of the two reactions and A{Z/A)'^ is 
the difference of the squared asymmetries of the sources 
in the two reactions. With A{Z/A)'^ and the tempera- 
ture determined as above, the free symmetry energy is 
extracted. 

From the free symmetry energy derived in this way 
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FIG. 1: (Color online) Free and internal symmetry energy 
as a function of the surface velocity. Experimental results 
are compared with results of theoretical calculations neglect- 
ing cluster formation (RMF) and including cluster formation 
(QS). 



from the measured yields, the internal symmetry energy 
can be calculated if the symmetry entropy is known. The 
values of the symmetry entropy S^J^ for given param- 
eters of temperature and density within the NSE model 
are shown in Tab. |T1 column 5. They are calculated with 
the equivalent expression of Eq. (1) as the difference be- 
tween the entropies of pure proton or neutron and sym- 
metric nuclear matter. In contrast to the mixing entropy 
that leads to a larger entropy for uncorrelated symmet- 
ric matter in comparison with pure neutron matter, the 
formation of correlations, in particular clusters, will re- 
duce the entropy in symmetric matter, see also Fig. 9 
of Ref . Q . For parameter values for which the yields of 
free nucleons in symmetric matter are small, the sym- 
metry entropy may become positive, as seen in Tab. |l] 
for low temperatures. The fraction of nucleons bound in 
clusters can decrease, e.g. due to increasing temperature 
or the dissolution of bound states at high densities due 
to the Pauli blocking. Then, the symmetric matter re- 
covers its larger entropy so that the symmetry entropy 
becomes negative, as seen in Tab. |T] also in the QS and 
self-consistent (sc, see below) calculations. 

The results obtained in this way for the internal sym- 
metry energy i^sym = ^svm + TS^, ^ a re shown in Tab. HI 
column 6. We note that in Ref. [20[ the symmetry en- 
tropy was estimated using results of the virial expansion 
of Ref. [3l leading to different internal symmetry ener- 
gies. However, this approximation is unreliable at the 
densities considered here. 

In Tab. HI we also give results of the QS model Q for 
the free and internal symmetry energies (columns 7 and 
8) at given T and n. In Fig. [1] (left panel) the experimen- 
tally obtained free symmetry energy is compared to the 
results of the RMF calculation without clusters and the 
QS model with clusters @. There are large discrepancies 
between the measured values and the results of calcula- 
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FIG. 2: (Color online) Comparisons of the scaled internal 
symmetry energy Esymin) / Esyni{no) as a function of the 
scaled total density n/no for different approaches and the 
experiment. Left panel: The symmetry energies for the com- 
monly used MDI parametrization of Chen et al. 23] for T — 
and different asy-stiffnesses, controlled by the parameter x 
(dotted, dot-dashed and dashed (black) lines); for the QS 
model including light clusters for temperature T — 1 MeV 
(solid (green) line), and for the RMF model at T = includ- 
ing heavy clusters (long-dashed (orange) line). Right panel: 
The internal scaled symmetry energy in an expanded low- 
density region. Shown are again the MDI curves and the QS 
results for T — 1,4, and 8 MeV compared to the experimental 
data with the NSE entropy (solid circles) and the results of 
the self-consistent calculation (open circles) from Tab. [T] 



tions in the mean-field approximation when cluster for- 
mation is neglected. On the other hand, the QS model 
results correspond nicely to the experimental data. In 
the right panel of Fig. [1] we compare the internal symme- 
try energy derived from the experimental data with the 
RMF and QS results. Again, it is clearly seen that the 
quasiparticle mean-field approach (RMF without clus- 
ters) disagrees strongly with the experimentally deduced 
symmetry energy while the QS approach gives a rather 
good agreement with the experimental data. 

In Fig. [2] we present results for different approaches 
to extracting the internal symmetry energy and compare 
with the experimental values. In the left panel of the fig- 
ure we show theoretical results for T at or close to zero. 
A widely used momentum-dependent parametrization of 
the symmetry energy (MDI) at temperature T = MeV 
was given in Refs. [J, |23[ and is shown for different as- 
sumed values of the stiffness parameter x. For these 
parametrizations the symmetry energy vanishes in the 
low-density limit. We compare this to the QS result at 
T = 1 MeV (at lower T crystallization or Bose condensa- 
tion may occur as discussed above). In this approach the 
symmetry energy is finite at low density. The T = 1 MeV 
curve will also approach zero at extremely low densities of 
the order of 10~^ fm~'^ because the temperature is finite. 
The RMF, T = curve is discussed below. Also note 
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that the underlying RMF model for the quasiparticle de- particle descriptions. 



scription with uq = 0.149 fm^^^ Esym{no) = 32.73 MeV 
gives a reasonable behavior at high density similar to 
the MDI. x^Q paramctrization. We thus see that our 
approach successfully interpolates between the cluster- 
ing phenomena at low density and a realistic description 
around normal density. 

In the right panel of Fig. [2] we compare to the ex- 
perimental results, full (red) circles (Tab. HI col. 6) 
in an expanded low-density region. Besides the MDI 
paramctrization we show the QS results [1| for T=l, 4, 
and 8 MeV, which are in the range of the temperatures in 
the experiment. The QS results including cluster forma- 
tion agree well with the experimental data points, as seen 
in detail in Fig. [TJ We conclude that medium-dependent 
cluster formation has to be considered in theoretical mod- 
els to obtain the low-density dependence of the symmetry 
energy that is observed in experiments. 

The temperatures and densities of columns 2 and 3 in 
Tab. H] will be modified if medium effects on the light 
clusters are taken into account We have carried 

out a self-consistent determination of the temperatures 
T^'^ and densities n^'^ taking into account the medium- 
dependent quasiparticle energies as specified in Ref. [12 1 



(columns 9 and 10 of Tab.[T|. Compared to the Albergo 
method results the temperatures T^'^ are about 10 % 
lower. Significantly higher values are obtained for the in- 
ferred densities rf^ which are more sensitive to the inclu- 
sion of medium effects. We have also calculated the free 
and internal symmetry energies corresponding to these 
self-consistent values of T^'^ and n^'^ according to Ref. Q 
(columns 11 and 12 of Tab. |T|. These results are also 
shown in the right panel of Fig. 2 as open (purple) cir- 
cles. The resultant internal symmetry energies are 15 to 
20 % higher than the QS model values for T and n given 
in columns 2 and 3 in Tab. IH 

We have restricted our present work to that region of 
the phase diagram where heavier clusters with A > 4 are 
not relevant. The generalization of the given approach 
to account for clusters of arbitrary size would lead to 
an improvement in the low-density, low-temperature re- 
gion when nuclear statistical equilibrium is assumed. Al- 
ternatively, one can introduce the formation of heavier 
nuclei in the presence of a nucleon and cluster gas, cf. 

Rcfs. di,!!^. 

The simplest approach to model the formation of heavy 
clusters is to perform inhomogcneous mean-field calcu- 
lations in the Thomas-Fermi approximation assuming 
spherical Wigner-Seitz cells. In Fig. [5] (left panel) pre- 
liminary results for the zero-temperature symmetry en- 
ergy of such a calculation is shown by the long-dashcd 
line using the same RMF paramctrization as for the QS 
approach introduced above; for details see Ref. The 
symmetry energy in this model approaches a finite value 
at zero density in contrast to the behavior of the MDI 
parametrizations and conventional single- nucleon quasi- 



In conclusion, we have shown that a quantum- 
statistical model of nuclear matter, that includes the 
formation of clusters at densities below nuclear satura- 
tion, describes quite well the low-density symmetry en- 
ergy which was extracted from the analysis of heavy-ion 
collisions. Within such a theoretical approach the com- 
position and the thermodynamic quantities of nuclear 
matter can be modeled in a large region of densities, 
temperatures and asymmetries that are required, e.g., 
in supernova simulations. 
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